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A METHOD FOR ENHANCING THE EDGE CONTRAST OF A DIGITAL 
IMAGE INDEPENDENTLY FROM THE TEXTURE 

CROSS-REFERENCE TO RELATED APPLICATION 

5 Reference is made to commonly assigned copending application 

Serial No. 09/324,239, entitled "A Method for Enhancing the Edge Contrast of a 
Digital Image" and filed June 2, 1999 in the names of Andrew C. Gallagher and 
Edward B. Gindele, which is assigned to the assignee of this application. 

1 0 FIELD OF THE INVENTION 

The invention relates generally to the field of digital image 
processing and, more particularly, to a method for enhancing the edge contrast and 
texture of a digital image. 

1 5 BACKGROUND OF THE INVENTION 

Traditional methods of increasing the apparent sharpness of a 
digital image, such as the technique of unsharp masking, often produce unwanted 
artifacts at large transition edges in the image. For example, unsharp masking is 
often described by the equation: 

20 Sproc - Sorg+B(Sorg-Sus) 

where Sproc represents the processed image signal in which the high frequency 
components have been amplified, Sorg represents the original image signal, Sus 
represents the unsharp image signal, typically a smoothed image signal obtained 
by filtering the original image, and B represents the high frequency emphasis 

25 coefficient. 

The unsharp masking operation may be modeled as a linear system. 
Thus, the magnitude of any frequency in Sproc is directly dependent upon the 
magnitude of that frequency in the Sorg image signal. As a consequence of this 
superposition principle, large edges in the Sorg image signal will often display a 
30 ringing artifact in the Sproc signal when the desired level of high frequency 

enhancement has been performed in other areas of the Sproc signal. This ringing 
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artifact appears as a light or dark outline around the large edge, and may be 
visually objectionable. 

Many non-linear filters based on local statistics exist for the 
purposes of noise reduction, sharpening, and contrast adjustment. For example, the 
5 median filter is well known in the art. In this filter, typically implemented for noise 
reduction, each pixel is replaced with the median value of some surrounding 
neighborhood. This filtering process is generally very successful at removing 
impulse noise; however, the processed image appears slightly less sharp than the 
original. 

1 0 Another example of a non-linear filter based on local statistics is 

local histogram equalization, referred to as adaptive histogram modification by 
William Pratt on pages 278 - 284 of the book Digital Image Processing, Second 
Edition , John Wiley & Sons, 1991. With this filter, the values of pixels are 
modified by the cumulative histogram of a local window. This technique 

1 5 effectively adjusts the contrast of each region of a digital image, effectively 
increasing the local contrast in some regions of the image, and decreasing the 
contrast in other regions. This technique does not intend to increase the apparent 
sharpness of any given region. Also, this technique does not ensure that the typical 
artifacts of ringing will not occur. 

20 There exist many algorithms for sharpening the appearance of 

images without generating artifacts or enhancing the notability of noise. In U.S. 
Patent 4,571,635, Mahmoodi and Nelson teach the method of deriving an emphasis 
coefficient B that is used to scale the high frequency information of the digital 
image depending upon the standard deviation of image pixel values in a local 

25 neighborhood. In addition, in U.S. Patent 5,081,692, Kwon and Liang teach that 
the emphasis coefficient B is based upon a center weighted variance calculation. 
However, neither use separate strategies for texture and edge regions. 

In U.S. Patent 4,794,531, Morishita et al teaches a method of 
generating the unsharp image with a filter whose weights on neighboring pixels are 

30 based upon the absolute difference between the central pixel and the neighboring 
pixel. Morishita claims that this method effectively reduces artifacts seen at the 
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edges of a sharpened image (as compared with traditional unsharp masking). 
However, the method of Morishita does not offer explicit control over the amount 
of edge reshaping. 

Thus, there exists a need for an alternative method of manipulating 
5 a digital image in order to generate an image signal that appears to be sharper, or 
more in focus, while minimizing the ringing artifact that is evident with the 
unsharp masking technique. 

SUMMARY OF THE INVENTION 

10 It is an object of this invention to allow for independent control of 

the enhancement applied to detail, large edges, and noisy regions. 

The present invention is directed to overcoming one or more of the 
problems set forth above. Briefly summarized, according to one aspect of the 
present invention, the invention resides in a method for enhancing a digital image 

1 5 channel comprised of a plurality of image pixels, particularly where the digital 
image channel is split into first and second signals, such as pedestal and texture 
signals, that substantially comprise the digital image channel. The method 
involves utilizing a predetermined tone scale conversion to enhance the digital 
image channel. Initially, image values are provided from the first signal, e.g., the 

20 pedestal signal, corresponding to image pixels from a region of the image. Then, 
a statistical characteristic of the image pixels in the region is identified, and the 
predetermined tone scale conversion is normalized for the statistical characteristic 
in order to generate a normalized tone scale conversion. The normalized tone 
scale conversion is then performed on a central pixel of the region in order to 

25 generate a first signal with enhanced image values; and the first signal with 
enhanced image values is combined with the second signal to generate an 
enhanced digital image channel. 

The present invention has the advantage of controlling the value of 
the central pixel in the pedestal signal in accordance with a statistical 

30 characteristic of the region, e.g., driving the value of the central pixel either 

toward the local maximum or the local minimum of the region, except in the case 
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where the central pixel is substantially midway between the local maximum and 
the local minimum. Consequently, edge transitions in the pedestal signal occur 
over a narrower range of pixels than in the input image, thus when combined with 
the texture signal generating an image signal that appears to be sharper, or more in 
5 focus, than the original image. Moreover, since the output of the tone scale 

conversion is modified by the statistical characteristic, e.g., bounded by the local 
maximum and the local minimum of the region, systematic overshoot and 
undershoot at an edge boundary is diminished and the ringing artifact is not as 
noticeable. 

1 0 These and other aspects, objects, features and advantages of the 

present invention will be more clearly understood and appreciated from a review 
of the following detailed description of the preferred embodiments and appended 
claims, and by reference to the accompanying drawings. 

1 5 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram illustrating an overview of the present 

invention. 

Fig. 2 is a block diagram of a preferred embodiment of the 
invention, showing the preference sharpener in Fig. 1 in more detail. 
20 Fig. 3 is a block diagram of a first embodiment of the sharpening 

processor shown in Fig. 2. 

Fig. 4 is a block diagram of a second embodiment of the sharpening 
processor shown in Fig. 2. 

Fig. 5 is a block diagram of a first embodiment of the pedestal 
25 splitter shown in Fig. 3. 

Fig. 5 A is a block diagram of a second embodiment of the pedestal 
splitter shown in Fig. 3. 

Fig. 6 is a block diagram of the avoidance signal generator shown in 

Fig. 5. 

30 Fig. 7 is n example of an artifact avoidance function applied by the 

avoidance signal generator shown in Fig. 6. 
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Fig. 8 is a block diagram of the pedestal modifier shown in Fig. 3. 

Fig. 9 is an example of a portion of an image showing a central 
pixel and the associated local region of the image. 

Fig. 10 is an example of several tone scale functions. 
5 Fig. 1 1 is a block diagram of the texture modifier shown in Fig. 3. 

Fig. 12 is a block diagram of a first embodiment of the local SNR 
estimator shown in Fig. 1 1 . 

Fig. 13 is a block diagram of a second embodiment of the local 
SNR estimator shown in Fig. 1 1 . 
1 0 Fig. 14 is a plot of a function used to generate a noise table. 

Fig. 15 is a plot of a look up table (LUT) used by the gain 
determiner shown in Fig. 1 1 . 

DETAILED DESCRIPTION OF THE INVENTION 

15 In the following description, the present invention will be described 

as a method implemented as a software program. Those skilled in the art will 
readily recognize that the equivalent of such software may also be constructed in 
hardware. Because image enhancement algorithms and methods are well known, 
the present description will be directed in particular to algorithm and method steps 

20 forming part of, or cooperating more directly with, the method in accordance with 
the present invention. Other parts of such algorithms and methods, and hardware 
and/or software for producing and otherwise processing the image signals, not 
specifically shown or described herein may be selected from such materials, 
components and elements known in the art. Given the description as set forth in 

25 the following specification, all software implementation thereof is conventional 
and within the ordinary skill in such arts. 

It is instructive to note that the present invention utilizes a digital 
image which is typically a two-dimensional array of red, green, and blue pixel 
values or of single monochrome pixel values corresponding to light intensities. In 

30 addition, the preferred embodiment is described with reference to an image of 

1024 rows of pixels and 1536 lines of pixels, although those skilled in the art will 



-6- 



recognize that digital images of different resolutions and dimensions may be used 
with equal, or at least acceptable, success. With regard to matters of nomenclature, 
the value of a pixel of a digital image located at coordinates (x,y), referring to the 
x th row and the y th column of the digital image, shall herein comprise a triad of 
5 values [r(x,y), g(x,y), b(x,y)] respectively referring to the values of the red, green 
and blue digital image channels at location (x,y). In this regard, a digital image 
may be considered as comprising a certain number of digital image channels. In 
the case of a digital image comprising red, green and blue two-dimensional arrays, 
the image comprises three channels, namely, red, green and blue channels. 

1 0 Additionally, a luminance channel n may be formed from the color signals. The 
value of a pixel of a digital image channel n located at coordinates (x,y), referring 
to the x th row and the/ A column of the digital image channel, shall herein be a 
single value referred to as n(x,y). 

Referring to the overview of the invention shown in Fig. 1, a digital 

1 5 image is input to a preference sharpener 2, which sharpens the input digital image 
for the preference of an operator. The purpose of the preference sharpener 2 is to 
enhance the detail present in the digital image without enhancing noise or creating 
artifacts. This purpose is accomplished by a process that decomposes a digital 
image channel into a signal corresponding to the image detail, a signal 

20 corresponding to the image edges, and a signal corresponding to the image signal 
to noise ratio (SNR), as will be described in detail hereinbelow. This 
decomposition and creation of separate signals allows for independent control of 
the enhancement applied to detail, large edges, and noisy regions. The output of 
the preference sharpener 2 comprises a digital image that, in accordance with the 

25 invention, appears sharper and more natural than can be achieved with traditional 
means of sharpening a digital image. 

Fig. 2 shows a block diagram of the preference sharpener 2. The 
digital color image and noise information (in the form of a noise table 6) are input 
to the preference sharpener 2. The noise information is input in the form of a 

30 luminance table. The digital image is input to a luminance/chrominance converter 
10 in order to generate a luminance digital image channel n(x,y) and two 
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chrominance digital image channels g(x,y) and i(x,y). The luminance digital image 
channel n(x,y) and the luminance noise table are input over respective lines 11a 
and 1 lb to a sharpening processor 20 for enhancing the digital image channel 
n(x,y) based on knowledge of the system noise. The transformation performed in 
5 the converter 10 from the digital color image, typically in an RGB space 
comprising separate red, green and blue image channels, to a luminance- 
chrominance color space is done by means of a color space matrix transformation, 
which results in the luminance digital image channel n(x,y) and the two 
chrominance digital image channels g(x 5 y) and i(x,y),as is well known in the art. In 
1 0 the preferred embodiment, the matrix transformation to the luminance- 
chrominance space may be described with the following equation: 
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Equation Al provides a matrix rotation from RGB space into a luminance- 
1 5 chrominance space in order to isolate a single channel upon which the tone scale 
function (to be described) can operate. For example, a film scanner records (for 
each pixel) the red, green and blue printing density at that point. These three 
values can be thought of as a pixel's location in three-dimensional space. The 
converter 10 performs an axis rotation as shown in equation Al, which provides 
20 for a neutral axis, upon which R = G = B, and two color difference axes, green- 
magenta and illuminant. 

The noise table input to the sharpening processor 20 provides the 
relationship between the signal intensity level i and the expected amount of noise 
a N (0 for that intensity level. In the preferred embodiment, as will be described in 
25 more detail, the noise table is a two column table, wherein the first column 
represents an intensity level / and the second column represents the expected 
standard deviation of noise a n (z) for that intensity level. 

The chrominance digital image channels g(x,y) and i(x,y) are input 
over a line 1 lc to a chrominance processor 40 and may be adjusted as desired. For 
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instance, the action of the chrominance processor 40 may be to scale the 
chrominance image channels by a constant greater than 1.0 in order to increase 
their apparent saturation. The operation of the chrominance processor is not of 
particular relevance to the present invention, and consequently it will not be 
5 discussed further, except to note that, in the present embodiment, the output of the 
chrominance processor 40 is maintained identical to its input. 

The digital image channels output from the sharpening processor 20 
and the digital image channels output from the chrominance processor 40 are input 
to an RGB converter 30 for conversion back to a digital image composed of a red , 

1 0 green , and blue digital image channels. This conversion is again accomplished 
with a matrix rotation (i.e., the inverse of the previous color rotation matrix 
performed by the converter 10). Inverting a 3 by 3 matrix is well known in the art 
and will not be further discussed. The output of the RGB converter 30 is a digital 
image that has been sharpened for the preference of the operator. 

15 As described in the preferred embodiment, the sharpening processor 

20 operates upon only a luminance digital image channel. However, as an 
alternative embodiment, the sharpening processor could be applied to each of the 
red, green, and blue digital image channels. In this case, the digital color image 
signals (which are RGB signals) would be directly applied to the sharpening 

20 processor 20. 

In Fig. 3, which illustrates a block diagram of the sharpening 
processor 20, the digital image channel n(x,y) is shown to be divided into two 
portions by a pedestal splitter 50, that is, into a pedestal signal and a texture signal. 
The preferred embodiment of the pedestal splitter 50, shown in Fig. 5, outputs a 

25 texture signal ntxt(x,y) and a pedestal signal n pe d(x,y). The texture signal n txt (x,y) 
primarily contains the detail of the image, as well as the image noise. The pedestal 
signal n P ed(x,y) is a conceptually smooth signal, except at regions of the image 
containing large occlusion boundaries, where the pedestal signal contains the 
transition edge. An alternative embodiment of the pedestal splitter 50, shown in 

30 Figure 5A, essentially comprises a linear FIR filter. In this alternative embodiment, 
the pedestal signal n pe d(x,y) is essentially equivalent to a lowpass signal and the 



-9- 



texture signal ntxt(x,y) is essentially equivalent to a highpass signal. Since other 
spatial frequency divisions are possible, those skilled in the art will understand that 
the present invention will yield beneficial results for any number of definitions of 
texture signals and pedestal signals. 
5 Referring again to Fig. 3, the digital image channel n(x,y), the 

texture signal n txt (x,y) and the noise table 6 are input to a texture modifier 70. The 
purpose of the texture modifier 70 is to enhance the magnitude of detail in the 
scene in a noise sensitive manner. An estimate of the local signal/noise ratio (SNR) 
is made using information about the expected standard deviation of noise for a 

1 0 given intensity level as supplied by the noise table. This estimate of the local SNR 
is used to set the boost factor relating to the level of enhancement given to a local 
level of the texture signal. This procedure will be described in detail hereinbelow. 
The output of the texture modifier 70 is an enhanced texture signal n' txt (x ? y). 
Furthermore, the pedestal signal n pe d(x,y) is input to a pedestal modifier 60 for the 

1 5 purpose of increasing edge contrast to make the edges appear to have greater 

clarity and sharpness. The purpose of the pedestal modifier 60 is to enhance the 
image edges without producing artifacts. The method employed by the pedestal 
modifier 60 will be described in detail hereinbelow. The output of the pedestal 
modifier is an enhanced pedestal signal n r pe d(x,y). 

20 The outputs of the pedestal modifier 60 and the texture modifier 70 

are input to an adder 80, which adds the enhanced texture signal n^t^y) and the 
enhanced pedestal signal np e d(x,y) together in order to produce the digital image 
channel output from the sharpening processor 20. The enhanced digital image 
channel n'(x,y) may be expressed as 

25 

n'(x,y) - n' txt (x,y) + n' pe d(x,y) 

Fig. 4 shows an alternative embodiment of the sharpening processor 
20, illustrating that minor perturbations to the arrangement of the components 
30 often have only a slight effect on the output of the sharpening processor 20. In this 
regard, the pedestal splitter 50 outputs the texture signal n txt (x ? y) and an avoidance 
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signal a(x,y). The avoidance signal a(x,y) is an intermediate signal calculated by 
the pedestal splitter 50 for determining the location of edges in the image. The 
derivation of this signal will be described below. The value of the avoidance signal 
ranges from 0.0 to 1.0. Image locations (x,y) having a value of the avoidance 
5 signal a(x,y) = 0.0 correspond with edge regions of the image; conversely, image 
locations (x,y) having a value of the avoidance signal a(x,y) > 0.0 correspond with 
regions of image detail or noise. As mentioned before, it is an object of this 
invention to allow for independent control of the enhancement applied to detail, 
large edges, and noisy regions. Thus, in the digital image channel output from the 
1 0 sharpening processor 20, it is desired that the edges (where a(x,y) = 0) be 

enhanced by the pedestal modifier 60, and that the regions of image detail or noise 
(corresponding to those locations where a(x,y) > 0.0) be enhanced by the texture 
modifier 70. 

To this end, the digital image channel n(x,y) is input to the pedestal 
1 5 modifier 60 and as before the texture signal ntxt(x,y) is input to the texture modifier 
70. The two resulting outputs from the pedestal modifier 60 and the texture 
modifier 70, the enhanced texture signal n't xt (x 5 y) and the enhanced pedestal signal 
n ! P ed(x,y), are input to an avoidance adder 81. The avoidance adder 81 requires 
three inputs: the two signals n't xt (x,y) and n' pe d(x,y) to be summed, and the artifact 
20 avoidance signal (a(x,y)). The two signals to be summed undergo a transformation 
such that one of the signals to be summed is multiplied by a(x,y), and the other is 
multiplied by (l-a(x,y)). The two transformed signals are then added. The signal 
input which is scaled by a(x,y) is known as the "a(x,y) input" of the avoidance 
adder 81 and the signal input which is scaled by (l-a(x,y)) is known as the "(1- 
25 a(x,y))" input of the avoidance adder 81 . In the case of the present embodiment of 
the sharpening processor 20, the signal output from the avoidance adder can be 
expressed as: 

n'(x,y) = a(x,y)ntxt(x,y) + (l-a(x,y))n ped (x 5 y) 



30 
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Ref erring to Fig. 5, the digital image channel n(x,y) input into the 
pedestal splitter 50 is split into a high pass signal nh P (x,y) and a low pass signal 
ni p (x,y) by a frequency splitter 94. Although there are many known techniques for 
generating highpass and lowpass signals, the frequency splitter is preferably 
5 implemented with a Gaussian filter having a standard deviation (sigma) of 0.9 
pixels. The preferred value of the standard deviation of the Gaussian filter may 
vary with image size. The value of 0.9 pixels for the standard deviation was 
derived by optimizing the present invention with 1024 by 1536 pixel size images. 
This Gaussian filter is a two dimensional circularly symmetric low-pass filter 
1 0 whose filter coefficients may be derived by the following formula which is well 
known in the art: 



g(i j) = l/(sigma sqrt(27t)) exp[-(i 2 +j 2 )/(2 sigma 2 )] 
Where g(i j) = the Gaussian filter coefficient at the (i, j) ^ pixel 
1 5 sigma = the standard deviation of the Gaussian filter (0.9) 

71 = the constant approximately 3.1415. 



In the preferred embodiment, i and j range from -3 to +3, inclusive, for a total of 
49 coefficients for the filter g(ij). Note that another technique, well known in the 
20 art, exists for separating a Gaussian filter into horizontal and vertical components 
for subsequent application for the purpose of computational cost reduction. In 
either case, the frequency splitter 94 applies the Gaussian filter g(i j) to the digital 
image channel n(x,y) by the process of convolution, which is well known in the 
art. The convolution may be expressed as: 

25 

ni p (x,y) = EE n(x-i,y-j)g(i j) 



with the summations occurring over all i and j. The signal ni p (x,y) resulting from 
this convolution process is the low-pass signal. The low-pass signal is output from 
30 the frequency splitter 94. In addition, the high-pass signal n hp (x,y) is output from 
the frequency splitter 94 after having been derived by the following relationship: 
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nh P (x,y) = n(x,y) - ni p (x,y) 

The low-pass signal ni p (x,y) is input into an avoidance signal 
5 generator 104 for forming the avoidance signal a(x,y), as will be described in detail 
below. A texture signal generator 90 receives both the high-pass signal n hp (x,y) 
and the avoidance signal a(x,y) and both signals are multiplied therein for 
producing the texture signal n txt (x,y). Thus, the texture signal n txt (x,y) may be 
expressed as: 

10 

ntxt (x,y) = a(x,y) * n hp (x,y). 

The texture signal n txt (x,y) computed by the texture signal generator 90 is then 

output by the pedestal splitter 50. In addition, as shown in broken line in Fig. 5, 
1 5 the avoidance signal a(x,y) may optionally be output by the pedestal splitter 50, in 

particular to provide an input to the avoidance adder 81 in the alternative 

embodiment of the sharpening processor 20 shown in Fig. 4. 

A pedestal generator 100 receives the luminance digital image 

channel n(x,y) and the texture signal ntxt (x,y) and subtracts the texture signal from 
20 the luminance signal, which produces the pedestal signal n pe a(x,y). Thus, the 

pedestal signal may be expressed as: 

n P ed (x,y) = n(x,y) - n tet (x,y). 

25 The pedestal signal n pe a(x,y) computed by the pedestal signal generator 100 is then 

output from the pedestal splitter 50. 

An alternative embodiment of the pedestal splitter 50 (identified as 

50A) is shown in Fig. 5A. The luminance digital image channel n(x,y) input into 

the alternative embodiment of the pedestal splitter 5 OA is split into a high pass 
30 signal and a low pass signal by a frequency splitter 94 preferably with the 

Gaussian filter as described above having a standard deviation (sigma) of 0.9 
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pixels. The preferred value of the standard deviation of the Gaussian filter may 
vary with image size. The value of 0.9 pixels for the standard deviation was 
derived by optimizing the present invention with 1024 by 1536 pixel size images. 
The frequency splitter 94 applies the Gaussian filter g(i j) to the digital image 
5 channel n(x,y) by the aforementioned process of convolution. 

The signal tii p (x,y) resulting from this convolution process is the 
low-pass signal output from the frequency splitter 94. In addition, a high-pass 
signal nh P (x,y) is output from the frequency splitter 94 after having been derived by 
the following relationship: 

10 

n hp (x,y) - n(x,y) - n ]p (x,y) 



In this alternative embodiment of the pedestal splitter, the lowpass 
signal ni p (x,y) is output from the pedestal splitter as the pedestal signal n pe d(x,y). 

1 5 (i.e., in this case, the pedestal signal n pe d(x,y) is set equal to the lowpass signal 
ni p (x,y) rather than being derived with a control signal). Likewise, the texture 
signal ntxt(x,y) is set equal to the highpass signal n hp (x,y). 

This alternative embodiment of the pedestal splitter 5 OA simply 
performs a digital filtering operation by convolution and outputs the highpass and 

20 lowpass signals as the texture and pedestal signals, respectively. This embodiment 
of the pedestal splitter 50A is simpler than the preferred embodiment 50 because 
no control signal is used or calculated. However, the quality of the digital image 
output from the preference sharpener 4 employing this alternative embodiment is 
inferior to the quality of the digital image output from the preference sharpener 4 

25 employing the preferred embodiment of the pedestal splitter 50. 

Those skilled in the art will recognize that there are many types of 
operations that may be used to generate two signals that essentially sum to 
generate the digital image channel. The preferred embodiment of the pedestal 
splitter and the alternative embodiment of the pedestal splitter are both examples 

30 of operations that may be used to separate a digital image channel into two signals 
that essentially sum to generate the digital image channel. 
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Referring to Fig. 6, there is illustrated a block diagram of the 
avoidance signal generator 104. In this regard, the low pass signal ni p (x,y) is input 
into a non-directional squared gradient calculator 106 for producing a non- 
directional gradient signal. This calculation is performed by first calculating the 
5 difference between the pixel and its upper vertical neighbor, as well as the 
difference between the pixel and its horizontal neighbor to the right. The non- 
directional squared gradient is the sum of the squares of these two differences. 
The non-directional squared gradient signal, ndg(x,y) may be expressed by the 
following equation: 

10 

ndg(x,y) = [nip(x,y)-ni p (x-l,y) f + [n lp (x,y)-m p (x,y+l) f 

The value of ndg(x,y) is then mapped by an artifact avoidance function applicator 
108 to produce an artifact avoidance signal a(x,y) as the output of the avoidance 
1 5 signal generator 104. 

Referring to Fig. 6, this mapping is performed by passing the non- 
directional squared gradient signal, ndg(x,y), through an artifact avoidance 
function av(y) formed in the following manner. The preferred embodiment of the 
present invention utilizes an artifact avoidance function defined by: 

20 

av(y) = (1/2) (1 + COSCTiCy-^/CC^Co)) 
for y > C 0 and y < C 1? 
25 av(y) = 0 

fory >=C 1 and 

av(y) = 1 

30 

for y <= C 0 where C Q and C 1 are numerical constants. 
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The preferred values for C 0 and C { will vary according to the range 
of the input data. The input data range in the images used in the development of 
this preferred embodiment ranged from 0 to 4095. In this case, the preferred value 
for C 0 is 996, and the preferred value for C x is 8400. As an alternate method for 
5 constructing the artifact avoidance function is described in the following equation: 

av(y) = l~sqrt((y-C 0 )/(C r C 0 )) 

for y >= C 0 and y <= C v 

10 

av(y) = 1 

for y < C 05 and 

15 av(y) = 0 

for y > C r Fig. 7 shows an example plot of the function av(y). 

Note that in the preferred embodiment, the value of Cj is related to 
the size sigma (in standard deviations) of the Gaussian filter. The value of C x is 
20 inversely proportional to the square of sigma (in standard deviations) of the 
Gaussian filter. Preferably, C x may be determined by the relation: 

Cj = 6804/(sigma*sigma) 

25 Also, C 0 is related to C x by the equation: 

C 0 = 0.127 C r 18 

The artifact avoidance signal, a(x,y), is generated by the artifact 
30 avoidance function applicator 108. This is accomplished by applying the artifact 
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avoidance function av(y), to the non-directional gradient signal ndg(x,y), as shown 
in Figure 7. The mathematical result is described by the equation: 

a(x,y) = av ( ndg(x ? y) ). 

5 

The artifact avoidance function shown in Figure 7 is most efficiently implemented 
as a look-up-table (LUT). 

It is helpful to understand that the avoidance signal a(x,y) is an 
example of a control signal created from a spatially filtered version of the digital 
1 0 image channel. To begin with, a generalized linear spatial filter is described by the 
equation 

c(x,y) - E£d(x-i,y-j)g(i j) 

where the d(x-i ? y-j) values represent local pixel values surrounding the (x 3 y)th 
1 5 pixel and the g(i.j) values represent numerical coefficients of a digital filter that do 
not depend on the pixel values. A non-linear spatial filter herein described shall be 
defined as a spatial filter which cannot be described by the linear spatial filter 
equation. The output of applying a generalized control signal to an input signal 
bears a multiplicative relationship to the input signal. An example of applying a 
20 control signal is given by the equation: 

c(x,y) = h(x,y)d(x,y) 

where the d(x,y) values represent the input signal (x ? y)th pixel values and the 
25 h(x,y) values represent the (x,y)th pixel values of the control signal. The result of 
applying a control signal to an input signal falls into the general category of a non- 
linear spatial filter if the control signal is derived from a spatial filtered version of 
the input signal. The avoidance signal a(x,y) is an example of a control signal 
created from a spatially filtered version of the digital image channel. The texture 
30 signal ntx t (x 5 y) described in the preferred embodiment is an example of a non- 
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linear spatial filter produced with the application of a control signal and a digital 
filter applied to the image channel n(x,y). 

Referring now to Fig. 8, the pedestal modifier 60 shown therein is a 
sharpening filter algorithm that is adaptive to the morphology of a digital image 
5 channel, as based on one or more image characteristics calculated over a localized 
region of the digital image channel (as isolated by a two-dimensional window). 
The filtering algorithm includes a stage for normalizing the tone scale conversion 
with a scaling function for the image statistic(s) over the localized region of a 
digital image channel and then using the scaled tone scale to reshape the structure 

10 of edges in the digital image channel. More specifically, the shape of the local 
tone scale function is ordinarily selected before processing the image, but the 
exact scaling and translation of the local tone scale operation is determined from 
the statistics in any localized region. Accordingly, after the region is isolated, 
implementation of the algorithm includes identifying a desired statistical 

1 5 characteristic from the local region, performing a local tone scale conversion that 
is normalized in relation to the statistical characteristic, and mapping the central 
pixel of the local region through the tone scale conversion to generate an 
enhanced central pixel value. More detail of the pedestal modifier 60 is described 
in commonly assigned U.S. Patent Application Serial No. 09/324,239, entitled "A 

20 Method for Enhancing the Edge Contrast of a Digital Image" and filed June 2, 

1999 in the names of A.G. Gallagher and E.B. Gindele, and which is incorporated 
herein by reference. 

The pedestal modifier 60 may be implemented in several modes. 
More specifically, and particularly as shown in the block diagram of Figure 8, the 

25 pedestal modifier 60 is preferably practiced in a first mode by scaling the input 

pixel values to the local maximum and minimum for each region and then applying 
the scaled input values to a tone scale function with normalized coordinates; the 
output obtained by this conversion operation is then processed with an inverse of 
the scaling function (used to generate the input) to obtain the enhanced output 

30 value. Alternatively, in a second mode of implementing the pedestal modifier 60, a 
local tone scale is constructed having coordinates that are bounded by the local 
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maximum and minimum for each region, and then each central pixel value is 
mapped through the bounded coordinates of the tone scale function. This will 
directly produce the enhanced output value. In a third mode, only the input of the 
tone scale function is scaled with the inverse scaling function and the output pixel 
5 is processed with the inverse of the scaling function; in a fourth mode, only the 
output of the tone scale function is scaled and the input pixel is processed with the 
scaling function. In any mode, the shape of the tone scale function may remain 
fixed over the image; however, it may be desirable to vary the shape of the tone 
scale function dependent upon the statistical characteristics of the region. 

1 0 Furthermore, while it may be typical to apply the normalized tone scale operation 
to substantially all the pixels in the input digital image channel, it may be desirable 
to only apply the method to selected areas of the digital image channel. For details 
of the second, third and fourth modes, the reader is directed to the aforementioned 
Serial No. 09/324,239, which is incorporated herein by reference. These 

1 5 alternative modes provide a processed digital image channel identical to the 

processed digital image channel obtained by the pedestal modifier 60 disclosed in 
the preferred embodiment (i.e., the first mode, Fig. 8). 

Referring further to Fig. 8, there is illustrated a block diagram of a 
preferred embodiment for implementing the pedestal modifier 60 according to the 

20 first mode as set forth above. It also facilitates understanding to note that the 
preferred embodiment of the pedestal modifier 60 is described in reference to a 
relatively high resolution digital image channel having, e.g., dimensions of 1024 
pixels in height by 1536 pixels in width, although the invention is similarly 
operative with image channels having greater or lesser image resolutions. The 

25 choice of image resolution, however, may have a concomitant effect on the 

window size, i.e., while a 5 X 5 window has been found preferable for 1024 X 
1536 resolution, the choice of a higher resolution image sensor would ordinarily 
allow a larger window size, and vice versa for a lower resolution sensor. 

With reference also to Fig. 9, the value of each input central pixel 

30 110 of the digital representation of the digital image channel is first passed to a 
local statistics determiner 116 for calculating local statistics. The values of 
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surrounding pixels 114 included within a local region 1 12 are also passed to the 
local statistics determiner 116. The local region 112 and the input pixel 110 
constitute a spatial region of the digital image channel. In the preferred 
embodiment, the local region about each input pixel 1 10 is defined as a square 
5 neighborhood centered on the input pixel for which the local statistics are being 
determined. In the preferred embodiment, the local region is a square 
neighborhood with 5 pixels on each side. Figure 9 illustrates the local region 
implemented in the preferred embodiment. The pixels 114 included within the 
local region 1 12 are illustrated as the 24 pixels (indicated by hatching) surrounding 

1 0 the input pixel 110 passed to the local statistics determiner 116. Those skilled in 
the art will recognize that the present invention will be operable with a large 
variety of local region sizes and shapes, and the usage (in this description and as 
claimed) of the term "central" to describe the processed pixel merely refers to a 
pixel having a local surround of pixels, rather than a pixel at the centroid of its 

1 5 local region. 

In Figure 8, the value u of each central pixel 1 10 is input to the 
local statistics determiner 116 and a scaling function e valuator 120. In addition, 
the values of the pixels 1 14 in the local region are input to the local statistics 
determiner 1 16 in order to generate a statistical characteristic of the local region. 

20 A scaling function generator 118 uses the statistical characteristic to generate a 
scaling function s(x) that is applied to the scaling function evaluator 120 in order 
to normalize the value u of each input pixel 110 and thereby generate a first 
intermediate value L The statistical characteristic is also used by an inverse 
scaling function generator 122 to generate a function s -1 (x) that is the inverse of the 

25 scaling function s(x). The first intermediate value I is then applied to a tone scale 
applicator 124, which maps the intermediate value I through a tone scale function 
f(x) provided by a tone scale function generator 126 in order to generate second 
intermediate value I2. The output intermediate value I 2 from the tone scale 
applicator 124 is then applied to an inverse scaling function applicator 128, which 

30 generates a processed value u* in the system image space by using the inverse 
function s _1 (x) provided by the inverse scaling function generator 122. 
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More specifically, the local statistics determiner 116 calculates 
several characteristics that describe the image data contained in the local region 
112 and the input pixel 110. The local statistics could be any mathematical 
combination of the pixel values contained in the local region 112 and the input 
5 pixel 110. In the preferred embodiment, the local statistics determiner 116 

determines two local statistics, the maximum and minimum pixel values from the 
local region 112 and the input pixel 110. As an alternative, the local statistics 
determiner 116 could determine the mean of the n maximum pixel values, and the 
mean of the m minimum pixel values, where n and m are integers much smaller 

1 0 than the number of pixels contained in the local region 112, for example, n = m = 
3. As will be described, these local statistics are utilized in the generation of the 
scaling function s(x) and the inverse scaling function s'^x). 

The local statistics determined by the local statistics determiner 116 
are then passed to the scaling function generator 118 and the inverse scaling 

1 5 function generator 122. In the preferred embodiment, the scaling function may be 
represented with the following equation: 




where X represents the maximum pixel value from the local region 112 and N 
20 represents the minimum pixel value from the local region 112. The purpose of the 
scaling function is to modify the value of the central pixel 110 such that the 
modified pixel value may be used as an input to a tone scale function 130, as 
shown in Figure 10. In the preferred embodiment, this tone scale function 130 is 
applied as a look-up-table (LUT). In the preferred embodiment, the effect of the 
25 scaling function is to normalize (i.e., compress) the interval [N X] of input values 
to an interval [0 1], as shown in Fig. 10. 

In a similar manner, the inverse scaling function generator 122 
creates the inverse scaling function s _1 (u); the inverse scaling function s _1 (u) may 
be represented by the following equation: 



30 
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s - l (u) = u(X-N) + N 



The purpose of the inverse scaling function is to scale the output values obtained 
from the tone scale function 130 back to the interval [N X]. The scaling function 
5 and the inverse scaling function are indeed inverses, as it can easily be shown that 



function evaluator 120 in order to generate an intermediate value I = s(u), where I 
is the intermediate value, s(x) is the scaling function, and u is the value of the input 

1 0 pixel 110, The output of the scaling function evaluator 120 is the intermediate 
value I. The intermediate value I is then passed to the tone scale applicator 124. 
The purpose of the tone scale applicator 124 is to affect the contrast of the value of 
the input pixel 110 with respect to the local region 112. The tone scale applicator 
124 modifies the intermediate value I, thereby creating the second intermediate 

1 5 value I2. The tone scale applicator 124 modifies the intermediate value I with the 
tone scale function 130 (i.e., f(x)). Thus, the operation of the tone scale applicator 
may be expressed in the following equation: I 2 = f(I), where 1 2 is the second 
intermediate value, f(x) is the tone scale function, and I is the intermediate value 
output from the scaling function evaluator 120. 

20 The tone scale function f(x) is generated by the tone scale function 

generator 126. In the preferred embodiment, the tone scale function f(x) is the 
sigmoid function, formed by integrating a Gaussian probability distribution 
function, and is well known in the art. The sigmoid function is fully characterized 
by the standard deviation a of the Gaussian probability function that is integrated. 

25 The tone scale function f(x) may be expressed by the following formula: 



s" i (s(u)) = u. 



The value of the input pixel 1 10 is then passed to the scaling 




dx 
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In the preferred embodiment, the tone scale function is found by computing the 
summation of discrete samples of a Gaussian distribution, as shown in the 
following formula: 



-x 



u-a{2x-\) 0 2 

exp 
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where 2a is the number of samples in the LUT. In the preferred embodiment, a = 
50. 



Notice that the slope of the tone scale function f(x) evaluated at x = 
0.5 may be determined by the relationship: 



1 




where y is the slope of the sigmoid function evaluated at x = 0.5, and the 
approximate slope of the summation described above, and a is the standard 
deviation of the Gaussian probability function. When the value of y> 1.0, the 
operation of the present invention sharpens the digital image channel by reshaping 
the local structure of edges. When the value of y is < 1.0, the operation of the 
present invention tends to produce an image with a softer appearance than the 
input image. In the preferred embodiment y = 2.1 and a = 0.19. In the preferred 
embodiment, the tone scale function f(x) is independent of the local statistics 
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output from the local statistics determiner 116, and thus need be computed only 
one time per image channel, despite the number of pixels being processed. 

Figure 10 illustrates a variety of tone scale functions, one of which 
is shown as the tone scale function 130, and the associated standard deviations of 
5 the Gaussian probability function c. Those skilled in the art will recognize that a 
wide variety of functions in addition to the sigmoid function will accomplish the 
goal of the tone scale applicator 124. 

The second intermediate value I2 output from the tone scale 
applicator 124 is passed to the inverse scaling function applicator 128. In addition, 

1 0 the inverse scaling function s _1 (x) output from the inverse scaling function 
generator 122 is passed to the inverse scaling function applicator 128. The 
purpose of the inverse scaling function applicator 128 is to generate the processed 
pixel value u' from the second intermediate value I2. The inverse scaling function 
applicator 128 generates the enhanced pixel value u T by evaluating the inverse 

1 5 scaling function: 



u^s' 1 ^) 
Substituting for I 2 and I, 

u'^fMu)}) 



20 



where s -1 (x) represents the inverse scaling function, f{x} represents the tone scale 
function, s(x) represents the scaling function, u represents the value of the input 
25 pixel 110, and u r represents the value of the enhanced pixel. The digital image 
channel output from the pedestal modified 60 is an enhanced pedestal signal 
n ! P ed(x,y). 

Typically, the process for generating a processed pixel value u f from 
the pedestal modifier 60 is repeated for every pixel of the digital image channel in 
30 order to generate a processed pixel value corresponding to each pixel in the digital 
image channel. These processed pixel values, when considered as a collection, 
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form a processed digital image channel. Alternatively, the process for generating a 
processed pixel value u 1 may repeated for a subset of every pixel of the digital 
image channel in order to generate a processed pixel value corresponding to a 
subset of pixels in the digital image channel. In the latter case, the subset, when 
5 considered as a collection, forms a processed image for a portion of the channel. In 
general, the parameters of the present algorithm may be adjusted in order to 
produce an enhanced digital image channel, with the appearance of greater edge 
contrast and more clarity than the digital image channel. 

The pedestal modifier 60 requires definition of two parameters: the 

1 0 size of the local region (which may be described by a window size for the 

sharpening filter) and the shape of the tone scale conversion function (which may 
be described by the slope of the tone scale function). The algorithm tends to 
perform best when the window size is kept fairly small, e.g., no greater than say 
about 7X7 pixels for a 1024 X 1536 imager. Indeed, as the window size and 

1 5 slope of the local tone scale function increase, the image begins to have a "paint- 
by-number" artifact. Also, a large window tends to compress texture on either side 
of a large occlusion boundary type edge in an image because the maximum and the 
minimum tend to occur on opposite sides of the edge. This effect is reduced by 
keeping the window size relatively small. A window size of 5 X 5 or fewer pixels 

20 has been found generally preferable for producing satisfactory results in typical 
situations for a 1024 X 1536 resolution. As mentioned earlier, other suitable 
window sizes may be empirically determined without difficulty for other image 
resolutions; alternatively, a larger window may be acceptable in certain situations 
for a given resolution if the user is willing to accept increased artifacts in order to 

25 obtain the sharpening effect of the larger window. 

The pedestal modifier as herein described has the advantage of 
controlling the value of the central pixel in accordance with a statistical 
characteristic of the region, e.g., driving the value of the central pixel either toward 
the local maximum or the local minimum of the region, except in the case where 

30 the central pixel is substantially midway between the local maximum and the local 
minimum. Consequently, edge transitions occur over a narrower range of pixels 
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than in the input image, thus generating an image that appears to be sharper, or 
more in focus, than the original image. Moreover, since the output of the tone 
scale conversion is modified by the statistical characteristic, e.g., bounded by the 
local maximum and the local minimum of the region, systematic overshoot and 
5 undershoot at an edge boundary is diminished and the ringing artifact is not as 
noticeable. 

A detailed view of the texture modifier 70 can be seen in Fig. 1 1 . 
The texture modifier 70 inputs the digital image channel n(x,y), the texture signal 
ntxt(x,y) and the luminance noise table 6. The output of the texture modifier 70 is 

10 an enhanced texture signal n'txt(x,y). The purpose of the texture modifier 70 is to 
boost the texture signal by a factor that is related to an estimate of the local signal 
to noise (SNR) ratio. Thus, the portion of the texture signal coincident with regions 
of the digital image channel having high SNR will experience a greater level of 
boost as compared with the portion of the texture signal associated with low SNR 

1 5 regions of the digital image channel. For example, suppose that an image to be 
sharpened contains a large region of clear blue sky. The only modulation in this 
region is likely to be noise resulting from the imaging system. Therefore, it is not 
desirable to increase to amplitude of this modulation. On the other hand, 
modulation of a grassy field is a result of the modulation due to many blades of 

20 grass in addition to image system noise. In this case, it is preferable to boost the 
modulation since it can be attributed to actual modulation in the scene. In the 
present invention, the signal to noise ratio is estimated to act as a classifier that 
may be used to distinguish between, e.g., the two aforementioned types of regions 
in a scene. 

25 In this regard, the digital image channel n(x,y), the texture signal 

ntxt(x,y) and the luminance noise table 6 are input to a local SNR estimator 160. 
The purpose of the local SNR estimator 160 is to estimate the local signal to noise 
ratio which will later be used to determine the level of boost given to a specific 
pixel of the texture signal. In the preferred embodiment, the output of the local 

30 SNR estimator 160 is one estimate of the SNR for each pixel in the digital image 
channel, represented by the symbol SNR(x,y). However, it is also possible that the 
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SNR may be calculated for only a subset of the pixels' locations of the digital 
image channel. The local SNR estimator 160 will be described in greater detail 
hereinbelow. 

The output of the local SNR estimator 160, SNR(x,y), is passed into 
5 a gain determiner 134. The gain determiner 134 is meant to determine an 

individual level of boost B(x,y) to apply at each particular (x,y) location of the 
texture signal, given that an estimate of the local SNR is known. This is 
accomplished by the use of a look up table (LUT), for example the LUT shown as 
a plot in Fig. 15. In the preferred embodiment, this LUT is monotonically 
1 0 increasing. The LUT output (the gain) increases with the increase of the LUT input 
(the estimate of the SNR.) 

The equation for the LUT shown in Fig 1 5 is given as: 



B(x,y) = sf™ for SNR(x,y) <SNR min 

1 5 B(x,y) = sf min +(sf max -sf mi ^ 

for SNR mi n<SNR(x,y)<SNR max 
B(x,y) = sf max for SNR(x,y) > SNR max 



Where sfmax is > sf mm? and SNR max ^ SNR mm . In the preferred embodiment, sf max is 

20 3.5, sf min is 2.0, SNR max is 2.0 and SNR min is 0.4. 

The gain B(x,y) output by the gain determiner 134 and the texture 
signal ntxt(x,y) are then both passed to a texture booster 136. The purpose of the 
texture booster 136 is to multiply each value of the texture signal ntxt(x,y) by its 
associated gain. The output of the texture booster 136 is the enhanced texture 

25 signal n' txt (x,y), given by the equation: 



n' txt (x,y) = n tx t(x ? y) * B(x,y). 



30 



This enhanced texture signal n f txt (x ? y) output from the texture modifier 70 is then 
passed to the adder 80 in Fig. 3 (or to the avoidance adder 81 in the case of the 
alternative embodiment shown in Fig. 4.) Note that when B(x,y) < 1.0, the 
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magnitude of the enhanced signal is less than the magnitude of the original texture 
signal |n tx t(x ? y)[>|nW(x ? y)|. Thus, if the present invention is tuned such that in 
noisy regions, B(x,y) <1.0, the processed digital image channel will appear to have 
less noise than the digital image channel input to the preference sharpener 2. 
5 A detailed block diagram of the local SNR estimator 160 may be 

seen in Fig 12. The digital image channel n(x,y) is first input to a local variance 
computer 170. The purpose of the local variance computer 170 is to attain a local 
estimate of signal activity. In the preferred embodiment, the local variance 
computer 170 calculates the local variance a R 2 (x,y) over a window of size W by W 
1 0 pixels. In the preferred embodiment W = 7 pixels, although it has been 

determined that alternately sized windows perform the task with nearly equivalent 
results. 

The computation of the local variance of a set of values is well 
known in the art, and may be expressed as: 



15 



■2<l 



Where i and j both range from -3 to +3, inclusive, for a total of W 2 = 49 pixels 

included in the computation of the variance a R 2 (x,y). Calculation of the variance at 
20 edges involves mirroring, a technique well known in the art of image processing. 

The output of the local variance computer 170 is a local variance for each pixel in 

the digital image channel. 

A subtracter 172 inputs both the digital image channel n(x,y) and 

the texture signal n txt (x,y). The output of the subtractor 172 is the pedestal signal 
25 n pe d(x,y). The subtractor 172 determines the pedestal signal by subtracting the 

texture signal n txt (x,y) from the digital image channel n(x,y), according to the 

formula: 

n pe d(x,y) = n(x,y) - n m (x,y) 
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Note that this pedestal signal n pe d(x,y) is identical to the pedestal signal determined 
by the pedestal splitter 50. In a practical implementation of the present invention, 
the pedestal signal need only be computed once. 

The noise table 6 and the pedestal signal n pe d(x,y) are input an 
5 expected noise determiner 190. The noise table contains the relationship between 
expected noise and intensity i of the digital image channel n(x,y). In the 
preferred embodiment, the noise table is a two column table, a sample of which is 
shown hereinbelow in a two column list. The first column represents intensity level 
and the second column represents the expected standard deviation of noise for that 
1 0 intensity level. 



i 
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11.5 


100 


11.5 


450 


11.9 


550 


12.0 


650 


12.1 


750 


12.2 
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12.4 
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12.3 


1150 
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1250 


12.3 


1350 


12.4 


1450 


12.4 


1550 


12.2 


1650 


11.9 


1750 


11.5 
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11.4 
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11.6 


2050 


12.1 
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2150 


12.8 


2250 


13.8 


2350 


15.0 


2450 


16.4 


2550 


17.9 


2650 


19.9 



The noise table may also be specified in terms of a four column 
table, where the first column is the intensity i, and the second to fourth columns 
5 correspond to the standard deviations of noise expected in the red, green, and blue 
channels, respectively, of the digital image channel input to the preference 
sharpener 2. In this case, the expected standard deviation of noise for intensity 
level i of the luminance channel n(x,y) is dependent upon the coefficients of the 
color rotation matrix used to generate the luminance channel, as given by: 

10 

c n (i) = sqrt((a r (i)*rcoef) 2 + (a g (i)*gcoef) 2 + (a b (i)*bcoef) 2 )) 

where a n (i) is the expected standard deviation of noise at intensity i for the 
luminance digital image channel n(x,y). Likewise, ov(i), cr g (i), <Tb(i) are the 

1 5 expected standard deviations of noise of the respective red, green, and blue digital 
image channels at intensity i, and rcoef, gcoef, and bcoef are the coefficients of the 
color rotation matrix used to generate the luminance digital image channel n(x,y). 
As previously mentioned, in the preferred embodiment all three of these 
coefficients equal 1/3. Consequently, the luminance channel given in the preferred 

20 embodiment is given by: 

a n (i) = l/3*sqrt(a r (i) 2 + a g (i) 2 + a b (i) 2 ) 

Thus a noise table corresponding to the luminance digital image channel n(x,y) 
25 may be either supplied or calculated by the expected noise determiner 190. 
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The expected noise determiner 190 outputs an intensity dependent 
estimate of the standard deviation of noise <y^(x f y) for the intensity level of the 
pedestal signal n pe d(x,y). The expected noise determiner 190 either performs a table 
look-up, which relates the intensity of the pedestal signal to the expected standard 
5 deviation of noise o^(x,y) 9 or calculates the expected noise, as described above. 
The output a^(x,y) of the expected noise determiner 190 may be represented as: 

vnfcy) = <SN(n pe d(x,y)). 

1 0 If the value of n p ^(x f y) is not coincident with a value of i for which a n (z) exists, 
then linear interpolation is used to determine a value for a n (n pe dOtj>)). Linear 
interpolation is well known in the art and will not be further discussed. A plot of a 
noise table can be seen in Fig. 14. 

The procedure for generating a noise table such as the one 

1 5 illustrated in Fig. 14 can be found in U.S. Patent No. 5,641,596 issued June 24, 

1997 to Gray et al. or in U.S. Patent No. 5,923,775 issued July 13, 1999 to Snyder 
et al., both of which are incorporated herein by reference. 

An SNR calculator 1 80 inputs both the expected standard deviation 
of noise a^(x,y) output from the expected noise determiner 190, and the local 

20 variance a^(x 9 y) output from the local variance computer 170. The SNR calculator 
180 estimates a local signal to noise ratio SNR(x,>?) by computing the ratio of the 
signal variance (assuming that the local variance a R 2 (x,>>) is the sum {in an RMS 
sense} of the variance due to signal and the variance due to noise.) Thus, the 
output of the SNR calculator 180 is an estimate of the local SNR at each pixel 

25 location , according to the following formula: 



SNR(x,>>) = sqrt( (jR 2 (x,y) - a N 2 (x,>>) ) / a N (xj/). 
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Note that before attempting the square root operation, the value of a^fay) - 
<5i/(x,y) is checked to ensure that it is positive. If it is not, the value of SNR (x,y) is 
set to equal 0.0. 

Notice that many other estimates of the signal to noise ratio may be 
5 formulated without deviating from the spirit of the present invention. For example, 
an alternative embodiment of the local SNR estimator 160 is shown in Fig. 13. 
The purpose of this embodiment is to provide an estimate of the local SNR by a 
means that is computationally less expensive that by that shown in Fig. 12. In this 
alternative embodiment, the noise and signal deviation are both assumed to be 

1 0 Gaussian distributed, which allows for the use of mean absolute deviation (MAD) 
rather than standard deviation. The calculation of mean absolute deviation does not 
involve the computationally expensive square root and squaring operations. In this 
embodiment, the texture signal n m (x y y) and the digital image channel n(x y y) are 
again input to the subtracter 172 for the purpose of producing the pedestal signal 

1 5 n pe d(x,y), as previously described. The pedestal signal n pe d(x,y) will again be used 
as an input to the expected noise determiner 190 to determine the value of a^fay)* 
The texture signal n txt (x,y) is input to the local MAD calculator 200 
in order to determine the local mean absolute deviation in the texture signal. By 
assuming that the mean value of the texture signal n txt (x,y) is 0, the formula for 

20 computing the local MAD used by the local MAD calculator 200 is given as: 



MAD R (x,y) = £| nUx+i ,y+j) I /W 2 



where i and j preferably both range from -3 to +3, for a total of W 2 = 49 pixels 
25 included in the calculation of the mean absolute deviation. The local MAD 

calculator 200 outputs the value of the mean absolute deviation MAD R (x,y), which 
will later be used in calculations necessary for determining the estimate of the local 
SNR. 

The purpose of a noise table adjuster 210 is to compensate for the 
30 differences between calculating MADr (x,y) of the texture signal, and calculating 
the standard deviation GR(x y y) of the digital image channel. The noise table adjuster 
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210 modifies the luminance noise table for the purpose of using mean absolute 
deviation, rather than standard deviation. Assuming, as previously mentioned, that 
the luminance noise table contains the relationship between expected noise 
and intensity i, It is known in the art that the MAD of a Gaussian distribution is 
5 related to the standard deviation a of that distribution by the following 
relationship: 

MAD =sqrt(2/7i)a 

1 0 Thus, in order to translate the luminance noise table from the metric of standard 
deviation to the metric of MAD, each value of a n (i) must be scaled by about 0.8. 

In addition, there is the difference between this alternative 
embodiment and the preferred embodiment concerning which signal is used for 
computation of the local signal activity measure. In the preferred embodiment, the 

1 5 variance of the digital image channel is computed. In this alternative embodiment, 
the texture signal is used to calculate the value of MAD R (x,y), providing a 
computational benefit when the signal mean is assumed to be zero. Thus the 
relationship of the texture signal n m (x,y) to the digital image channel n(x,y) must 
be understood in order to adjust the luminance noise table to account for this. This 

20 relationship can be easily expressed only if the digital filtering process used to 
generate the texture signal is linear. For the sake of this step, the effects of the 
artifact avoidance signal are ignored, and the texture signal is assumed to be 
equivalent to the highpass signal. For a general (n x m) filter with coefficients 
h(ij) 9 wherein i = -(n-l)/2, -(«-3)/2 v .. 9 «-l/2 and j = -(m-l)/2, -(m-3)/2,...,/w-l/2, the 

25 variance of the filtered signal is related to the variance of the original signal by the 
following relationship: 

a 2 fs - <y 2 0S YZh 2 (ij) 



30 



where a 2 fs represents the variance of the filtered signal, a 2 0S represents the 
variance of the original signal, and the summation occurs over i and j. 
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Thus, in the case of the present invention, the filter h(ij) is assumed 
to be expressed as: 



h(ij) =S(ij)-g(ij) 

5 

where 8(ij) is the delta function which has a value of one at i = 0 and j = 0 and a 
value of zero everywhere else. This expression accounts for the previously 
described process of generating the texture signal (ignoring the artifact avoidance 
signal) by blurring the digital image channel with a Gaussian filter g(ij) 9 then 
1 0 subtracting the result from the digital image channel. Thus, considering this 

description, the noise table adjuster 210 must scale each value of the luminance 
noise table by a factor to account for the estimation of local signal activity on the 
texture signal rather than the digital image channel. This factor/ can be expressed 
as: 

15 

/= sqrt (YE h(ijf) 



where the summation occurs over i and j 9 and h(ij) is as previously described. The 
square root of the result is taken because the luminance noise table is in the units 
20 of standard deviation, rather than variance. 

Thus, the operation of the noise table adjuster 210 is to scale each 
entry a n (0 of the luminance noise table by a factor m equal to 



m —f* sqrt ( 21%) 

25 

in order to compensate for the calculation of MAD rather than standard deviation, 
and the calculation with the texture signal rather then the digital image channel. 

The modified noise table output from the noise table adjuster 210 is 
passed as an input to the expected noise determiner 190 for the purpose of 
30 attaining an estimate of the expected MAD from noise, MAD N (x,y). The expected 
noise determiner 190 uses the modified noise table to determine the value of 
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MAD N (x,y) for each value of n pe ^(x f y). The output of the expected noise 
determiner 190 is the value of MAX>N(x,y) at each location of the digital image 
channel. 

The SNR calculator 180 of the alternative embodiment shown in 
5 Fig. 13 acts essentially the same as the SNR calculator 180 of the preferred 
embodiment. The purpose of the SNR calculator 180 is to output a value of 
estimated SNR at each location of the digital image channel. To this end, the SNR 
calculator 180 calculates the following formula: 



1 0 SNR(x,y) = sqrt(MAD R (x,y) 2 - MAD N (x,y) 2 ) / MAD N (x,y) 

The estimate of the local SNR is output from the local SNR estimator 160 and 
input to the gain determiner 134, shown in Fig. 1 1 in order to allow for the 
determination of the appropriate boost factor for the texture signal at each location. 

15 In summary, the use of the pedestal splitter 50, decomposing the 

image into a signal containing the image detail and another signal primarily 
containing the large image edges, allows for independent control of the 
enhancement applied to edges and texture. Edges are reshaped by a morphological 
operation in order to increase edge contrast and avoid boundary artifacts. Texture 

20 in enhanced with regard to an estimate of the local signal to noise ratio. 

The present invention has been described with reference to a 
preferred embodiment. Changes may be made to the preferred embodiment 
without deviating from the scope of the present invention. For instance, in the 
preferred embodiment, a calculation of the local SNR is made for each pixel of the 

25 digital image channel. However, it is conceivable that the calculation of the SNR 
could be made for only every N th pixel in order to reduce computational cost, and 
replicated or interpolated. Calculation of the local SNR was described with 
reference to a local variance calculation on the digital image channel n(x,y). 
Additionally, an alternative embodiment described a method by which local SNR 

30 may be estimated based upon using a MAD calculation upon the texture signal 
n tx t(x,y). Those skilled in the art will recognize that the local SNR may also be 
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estimated by performing a variance calculation on the texture signal, or by 
performing a MAD calculation upon the digital image channel n(x,y). Such 
modifications to the preferred embodiment do not significantly deviate from the 
scope of the present invention. 
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PARTS LIST 
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170 
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180 
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190 


expected noise determiner 


200 


local MAD calculator 


210 


noise table adjuster 
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WHAT IS CLAIMED IS: 

1 . A method for enhancing a digital image channel comprised of a 
plurality of image pixels, said method comprising the steps of: 

(a) splitting the digital image channel into first and second 
signals that substantially comprise the digital image channel, said first and second 
signals each characterized by image values corresponding to the image pixels; 

(b) providing a predetermined tone scale conversion to enhance 
the digital image channel; 

(c) providing image values from the first signal corresponding 
to image pixels from a region of the image: 

(d) identifying a statistical characteristic of the image pixels in 

the region; 

(e) normalizing the predetermined tone scale conversion for the 
statistical characteristic in order to generate a normalized tone scale conversion; 

(f) performing the normalized tone scale conversion on a 
central pixel of the region in order to generate a first signal with enhanced image 
values; and 

(g) combining the first signal with enhanced image values and 
the second signal to generate an enhanced digital image channel. 

2. The method of claim 1 in which the first signal is a low pass 
signal and the second signal is a high pass signal. 

3. The method of claim 1 in which the first signal is a pedestal 
signal characteristic of edges in the digital image channel and the second signal is 
a texture signal characteristic of texture in the digital image channel. 

4. The method of claim 1 in which the second signal is modified 
before it is combined with the first signal. 
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5. The method of claim 1 in which the region of step (c) comprises 
providing image pixels from a two dimensional region distributed around the 
central pixel. 

6. The method of claim 5 in which the region of step (c) comprises 
a 5 X 5 or smaller array of image pixels. 

7. The method of claim 1 in which step (d) comprises determining 
a gradient value within the region. 

8. The method of claim 7 in which the gradient value is derived 
from a maximum and a minimum value for the region. 

9. The method of claim 8 in which the gradient value is computed 
from the difference between maximum and minimum image values of image pixels 
within the region. 

10. The method of claim 7 in which the gradient value is computed 
from the difference between average maximum and minimum values across 
respective pluralities of image pixels comprising portions of the region. 

11. The method of claim 1 in which step (b) includes the 
application of a predetermined tone scale function, step (e) comprises normalizing 
the value of the central image pixel input to the tone scale function with a scaling 
function based on the statistical characteristic, and step (f) comprises mapping the 
scaled value of the central pixel through the tone scale function to generate an 
intermediate output value. 

12. The method of claim 1 1 in which step (f) further comprises 
processing the intermediate value of the central image pixel output from the tone 
scale function with a function that is the inverse of the scaling function. 
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13. The method of claim 1 in which step (b) includes the 
application of a predetermined tone scale function and step (e) comprises 
normalizing the coordinates of the tone scale function based on the statistical 
characteristic. 

14. The method of claim 1 in which step (b) includes application of 
a predetermined tone scale function, and the method further comprises the step of 
adjusting a parameter of the tone scale function based on the statistical 
characteristic of the region. 

15. The method of claim 14 in which the step of adjusting a 
parameter of the tone scale function comprises varying the slope of the tone scale 
function. 

16. The method of claim 1 in which the steps (c) - (f) are repeated 
for consecutive overlapping regions such that contiguous central pixels are 
processed through the tone scale conversion. 

17. The method of claim 16 wherein steps (e) - (f) include 
application of a predetermined tone scale function to all central image pixels 
processed through the normalized tone scale conversion. 

18. The method of claim 1 in which step (b) includes the 
application of a predetermined tone scale function having input and output 
coordinates, step (e) comprises normalizing the input coordinate of the tone scale 
function with a scaling function based on the statistical characteristic and 
normalizing the value of the central pixel output from the tone scale function based 
on the inverse of the scaling function, and then step (f) comprises mapping the 
value of the central pixel through the tone scale function to generate a scaled 
output pixel. 
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19. The method of claim 1 in which step (b) includes the 
application of a predetermined tone scale function having input and output 
coordinates, step (e) comprises normalizing the value of the central pixel input to 
the tone scale function with a scaling function based on the statistical 
characteristic and normalizing the output coordinate of the tone scale function 
based on the inverse of the scaling function, and then step (f) comprises mapping 
the scaled value of the central pixel through the tone scale function to generate a 
processed output pixel. 

20. A method for enhancing a digital image processed through a 
digital image channel comprised of a plurality of image pixels, the method 
comprising the steps of: 

(a) splitting the digital image into first and second signals that 
substantially comprise the digital image channel, said first and second signals each 
characterized by image pixel values corresponding to the image pixels; 

(b) receiving an image pixel value of the first signal from a spatial 
region of the digital image; 

(c) receiving a tone scale function; 

(d) calculating a statistical parameter based on the values of other 
pixels of the first signal in the spatial region of the digital image; 

(e) generating a scaling function based upon the statistical 

parameter; 

(f) generating a first intermediate value from the image pixel value 
and the scaling function; 

(g) generating a second intermediate value from the first 
intermediate value and the tone scale function; 

(h) generating an enhanced pixel value from an inverse of the 
scaling function and the second intermediate value; 

(i) combining the first signal with the enhanced pixel value and the 
second signal to generate an enhanced digital image channel. 
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21 . The method of claim 20 in which the first signal is a low pass 
signal and the second signal is a high pass signal. 

22. The method of claim 20 in which the first signal is a pedestal 
signal characteristic of edges in the digital image channel and the second signal is 
a texture signal characteristic of texture in the digital image channel. 

23. The method of claim 20 in which the second signal is modified 
before it is combined with the first signal. 

24. The method of claim 20 further comprising the step of 
repeating steps (b) and (d) - (i) in order to generate an enhanced pixel for each 
image pixel in the digital image channel. 

25. The method of claim 20 wherein the statistical parameter of 
step (d) is the maximum or minimum of the pixel values from the spatial region of 
the digital image channel. 

26. The method of claim 25 wherein the statistical parameter is 
used to normalize the tone scale function thereby generating a normalized tone 
scale function. 

27. The method of claim 26 wherein the statistical parameter is 
used to normalize the original pixel value creating a normalized pixel value. 

28. The method of claim 27 wherein the normalized pixel value is 
modified by the normalized tone scale function. 



29. The method of claim 27 wherein the tone scale function is a 
predetermined sigmoid function. 
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30. The method of claim 29 wherein the tone scale function is a 
predetermined sigmoid function used for each pixel in the digital image. 

3 1 . The method of claim 29 wherein the tone scale function is a 
predetermined sigmoid function that is normalized to create a normalized sigmoid 
function. 



32. The method of claim 3 1 wherein the normalized pixel value is 
modified by the normalized sigmoid function. 

33. The method of claim 29 wherein the step of using a sigmoid 
function is further comprised of determining the standard deviation of the sigmoid 
function based upon the statistical parameters. 

34. A computer program product for enhancing a digital image 
channel comprised of a plurality of image pixels, said computer program product 
comprising: a computer readable storage medium having a computer program 
stored thereon for performing the steps of: 

(a) splitting the digital image channel into first and second 
signals that substantially comprise the digital image channel, said first and second 
signals each characterized by image values corresponding to the image pixels; 

(b) providing a predetermined tone scale conversion to enhance 
the digital image channel; 

(c) providing image values from the first signal corresponding 
to image pixels from a region of the image: 

(d) identifying a statistical characteristic of the image pixels in 

the region; 

(e) normalizing the predetermined tone scale conversion for the 
statistical characteristic in order to generate a normalized tone scale conversion; 
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(f) performing the normalized tone scale conversion on a 
central pixel of the region in order to generate a first signal with enhanced image 
values; and 

(g) combining the first signal with enhanced image values and 
the second signal to generate an enhanced digital image channel. 

35. The computer program product of claim 34 in which the first 
signal is a low pass signal and the second signal is a high pass signal. 

36. The computer program product of claim 34 in which the first 
signal is a pedestal signal characteristic of edges in the digital image channel and 
the second signal is a texture signal characteristic of texture in the digital image 
channel. 

37. The computer program product of claim 34 in which the second 
signal is modified before it is combined with the first signal. 

38. The computer program product of claim 34 in which the region 
of step (c) comprises providing image pixels from a two dimensional region 
distributed around the central pixel. 

39. The computer program product of claim 38 in which the region 
of step (c) comprises a 5 X 5 or smaller array of image pixels. 

40. The computer program product of claim 34 in which step (d) 
comprises determining a gradient value within the region. 

41 . The computer program product of claim 40 in which the 
gradient value is derived from a maximum and a minimum value for the region. 
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42. The computer program product of claim 41 in which the 
gradient value is computed from the difference between maximum and minimum 
image values of image pixels within the region. 

43. The computer program product of claim 40 in which the 
gradient value is computed from the difference between average maximum and 
minimum values across respective pluralities of image pixels comprising portions 
of the region. 

44. The computer program product of claim 34 in which step (b) 
includes the application of a predetermined tone scale function, step (e) comprises 
normalizing the value of the central image pixel input to the tone scale function 
with a scaling function based on the statistical characteristic, and step (f) comprises 
mapping the scaled value of the central pixel through the tone scale function to 
generate an intermediate output value. 

45. The computer program product of claim 44 in which step (f) 
further comprises processing the intermediate value of the central image pixel 
output from the tone scale function with a function that is the inverse of the scaling 
function. 

46. The computer program product of claim 34 in which step (b) 
includes the application of a predetermined tone scale function and step (e) 
comprises normalizing the coordinates of the tone scale function based on the 
statistical characteristic. 

47. The computer program product of claim 34 in which step (b) 
includes application of a predetermined tone scale function, and the method further 
comprises the step of adjusting a parameter of the tone scale function based on the 
statistical characteristic of the region. 
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48. The computer program product of claim 47 in which the step of 
adjusting a parameter of the tone scale function comprises varying the slope of the 
tone scale function. 

49. The computer program product of claim 34 in which the steps 
(c) - (f) are repeated for consecutive overlapping regions such that contiguous 
central pixels are processed through the tone scale conversion. 

50. The computer program product of claim 49 wherein steps (e) - 
(f) include application of a predetermined tone scale function to all central image 
pixels processed through the normalized tone scale conversion. 

51 . The computer program product of claim 34 in which step (b) 
includes the application of a predetermined tone scale function having input and 
output coordinates, step (e) comprises normalizing the input coordinate of the tone 
scale function with a scaling function based on the statistical characteristic and 
normalizing the value of the central pixel output from the tone scale function based 
on the inverse of the scaling function, and then step (f) comprises mapping the 
value of the central pixel through the tone scale function to generate a scaled 
output pixel. 

52. The computer program product of claim 34 in which step (b) 
includes the application of a predetermined tone scale function having input and 
output coordinates, step (e) comprises normalizing the value of the central pixel 
input to the tone scale function with a scaling function based on the statistical 
characteristic and normalizing the output coordinate of the tone scale function 
based on the inverse of the scaling function, and then step (f) comprises mapping 
the scaled value of the central pixel through the tone scale function to generate a 
processed output pixel. 
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ABSTRACT 

A method is described for enhancing a digital image channel, 
particularly where the digital image channel is split into pedestal and texture 
signals that substantially comprise the digital image channel. The method utilizes 
5 a predetermined tone scale conversion to enhance the digital image channel. 
Initially, image values are provided from the pedestal signal corresponding to 
image pixels from a region of the image. Then, a statistical characteristic of the 
image pixels in the region is identified, and the predetermined tone scale 
conversion is normalized for the statistical characteristic in order to generate a 
1 0 normalized tone scale conversion. The normalized tone scale conversion is then 
performed on a central pixel of the region in order to generate a pedestal signal 
with enhanced image values; and the pedestal signal with enhanced image values 
is combined with the texture signal to generate an enhanced digital image channel. 
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SECOND GIVEN NAME 


0 


RESIDENCES 
CITIZENSHIP 


CITY 


STATE OR FOREIGN COUNTRY 


COUNTRY OF CITIZENSHIP 


3 


BUSINESS 
ADDRESS 


BUSINESS ADDRESS 


CITY 


STATE & ZIP CODE (COUNTRY) 


2 


FULL NAME OF 
INVENTOR 


FAMILY NAME 


FIRST GIVEN NAME 


SECOND GIVEN NAME 


0 


RESIDENCE & 
CITIZENSHIP 


CITY 


STATE OR FOREIGN COUNTRY 


COUNTRY OF CITIZENSHIP 


4 


BUSINESS 
ADDRESS 


BUSINESS ADDRESS 


CITY 


STATE & ZIP CODE (COUNTRY) 


2 


FULL NAME OF 
INVENTOR 


FAMILY NAME 


FIRST GIVEN NAME 


SECOND GIVEN NAME 


0 


RESIDENCE & 
CITIZENSHIP 


CITY 


STATE OR FOREIGN COUNTRY 


COUNTRY OF CITIZENSHIP 


5 


BUSINESS 
ADDRESS 


BUSINESS ADDRESS 


CITY 


STATE & ZIP CODE (COUNTRY) 


2 


FULL NAME OF 
INVENTOR 


FAMILY NAME 


FIRST GIVEN NAME 


SECOND GIVEN NAME 


0 


RESIDENCE & 
CITIZENSHIP 


CITY 


STATE OR FOREIGN COUNTRY 


COUNTRY OF CITIZENSHIP 


6 


BUSINESS 
ADDRESS 


BUSINESS ADDRESS 


CITY 


STATE & ZIP CODE (COUNTRY) 



Direct Telephone Calls to: 

(name and telephone number) 



David M. Woods 
(716) 477-5256 
FAX: (716)477-4646 



I hereby declare that all statements made herein of my own knowledge are true and that all statements 
true; and further that these statements were made with the knowledge that willful false statements 
imprisonment, or both, under section 1001 of Title 18 of the United States Code, and that such willful 
application or any patent issuing thereon. 



made on information and belief are believed to be 
and the like so made are punishable by fine or 
false statements may jeopardize the validity of the 
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SIGNATURE OF INVENTOR 202 



SIGNATURE OF INVENTOR 203 



SIGNATURE OF INVENTOR 204 



SIGNATURE OF INVENTOR 205 



SIGNATURE OF INVENTOR 206 
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